The dynamical frustration of interlayer excitons delocalizing in 
bilayer quantum antiferromagnets 
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Abstract - Using the self-consistent Born approximation we study the delocalization of inter- 
layer excitons in the bilayer Heisenberg quantum antiferromagnet. Under realistic conditions we 
find that the coupling between the exciton motion and the spin system is strongly enhanced as 
compared to the case of a single carrier, to a degree that it mimics the confinement physics of 
carriers in Ising spin systems. We predict that the 'ladder spectrum' associated with this con- 
finement physics should be visible in the c-axis exciton spectra of insulating bilayer cuprates such 
as YBa2Cu306. Our discovery indicates that finite density systems of such excitons should show 
very rich physical behavior. 



Introduction. — The discovery of high Tc supercon- 
ductivity triggered a concerted theoretical effort aimed at 
understanding the physics of doped Mott insulators [l]|2]. 
Although much is still in the dark, the problem of an iso- 
lated carrier in the insulator is regarded as well understood 
[sjjs]. It turned out to be a remarkable affair, rooted in 
the quantum-physical conflict between the antiferromag- 
netism of the spin system and the delocalizing carrier. 
This conflict is at its extreme dealing with a classical Ising 
spin system, where a famous cartoon arises for the idea 
of confinement (Fig. [l^): the hopping causes a 'magnetic 
string' of overturned spins between the delocalizing charge 
and the spin left at the origin with an exchange energy in- 
creasing linearly in their separation. It was realized that 
the quantummechanical nature of the 5 = 1/2 Heisenberg 
spin system changes this picture drastically. The quan- 
tum spin-fluctuations repair efficiently this 'conflnement 
damage' in the spin background and one finds a 'spin- 
liquid polaron' as quasiparticle that propagates coherently 
through the lattice on a scale set by the exchange constant. 
This physics can be reliably addressed by parametrizing 
the spin system in terms of its linear spin waves (LSW), 
while the strong coupling between the spin waves and the 
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propagating hole is well described in terms of the self con- 
sistent Born approximation (SCBA). This turned out to 
be accurate to a degree that the photoemission results in 
insulating cuprates were quantitatively explained in this 
framework [o]. 

A related problem is the delocalization of a bound 
electron-hole pair (exciton, or more exactly the bound 
state of a double occupied and vacant site) through the 
antiferromagnetic background. It is easy to see that the 
propagation of an exciton in a single layer is barely affected 
by the antiferromagnetism since the combined motion of 
the electron and the hole neutralize the 'damage' in the 
spin system ^lOj. A problem of contemporary interest is 
the exciton formed in a bilayer, where the electron and the 
hole reside in the different layers. Using modern growth 
techniques it appears possible to engineer systems where 
such excitons are formed at a finite chemical potential, 
making it possible to address for instance exciton Bose 
condensation in the background of an antiferromagnetic 
spin system. 

An exciton in a Mott insulator is, unlike its counterpart 
in semiconductors, the bound state of the double occu- 
pied and vacant sites. We specialize to the case where 
this Coulomb interaction is large compared to the kinetic 
energy scale. This might be a realistic assumption in for 
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Fig. 1: a. A moving hole in a classical single layer Ising anti- 
ferromagnet causes a string of spin frustration, which leads to 
localization, b. Zero temperature phase diagram of the bilayer 
Heisenberg model as a function of interlayer coupling strength 
a = ^ on the horizontal axis. At a critical value ac a quan- 
tum phase transition exists from the antiferromagnetic to the 
singlet phase. The vertical axis shows the Neel order parame- 
ter signaling antiferromagnetism. (Adapted from [l2].) c. In 
a perfect Neel bilayer, as seen here from the side, the motion 
of an exciton causes a mismatch in the spin ordering. Here the 
exciton has moved one position to the right, changing places 
with the electrons residing on those sites 



example the cuprates. Due to Coulomb repulsion between 
electrons the lowest energy state is always the state where 
the empty and double occupied site are nearest neighbor 
along the interlayer rung. Here we report the discovery 
that under these realistic conditions such bilayer excitons 
couple extremely strongly through their quantum motion 
to the spin system. 

In fact, when the interlayer exchange coupling is small 
and the exciton hopping rate is large, one enters a regime 
that is similar to the confinement associated with the Ising 
spins, although the spin system is in the quantized Heisen- 
berg regime. This is illustrated by the exciton spectral 
function shown in Fig. [2] as computed with the LSW- 
SCBA method, showing the non-dispersive 'ladder spec- 
trum' which is a fingerprint of confinement. Fig. [l]^ de- 
picts a cartoon of the confinement mechanism: every time 
the exciton hops it creates two spin flips in the different 
layers that can only be repaired by quantum spin fluctu- 
ations driven by the interlayer exchange coupling. The 
rapid intralayer quantum spin fluctuations are now inef- 
fective, because the restoration of the anti-ferromagnetism 
requires quantum spin fluctuations that occur simultane- 
ously in the two layers with a probability that is strongly 
suppressed. 

This confinement effect can be studied directly in exper- 
iment by measuring the exciton spectrum in c-axis optical 
absorption of the YBa2Cu306 (YBCO) insulating bilayer 
system. Using realistic parameters we anticipate that this 
will look like Fig. [3j the main difference with Fig. [2]is that 
the exciton hopping rate is now of order of the exchange 
energy and in this adiabatic regime the spectral weight in 
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Fig. 2: Exciton spectral function for J = 0.2t and a — 0.2. 
On top of the incoherent bump a strong ladder spectrum has 
developed, signaling Ising confinement. The exact Ising ladder 
spectrum is shown in green dotted lines. The Ising peaks are 
very weakly dispersive, with bandwidth of order J. 



the ladder spectrum states is reduced. Our discovery is 
particularly significant in the context of the effort to cre- 
ate bilayer exciton systems at a finite density. Stripes and 
other complex ordering phenomena in the finite density 
electron systems are believed to be driven by the hole- 
spin system 'quantum frustration'. Since the excitons are 
bosons, this holds the promise that the same type phenom- 
ena can now be studied in a much more tractable bosonic 
setting. 

The exciton i — J model. — Let us now discuss our 
calculations. Our point of departure is the bilayer Heisen- 
berg model 



(1) 
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where Su is the spin operator on site i and layer / = 1,2. 
This model is known to describe the physics of half-filled 
bilayer Mott insulators, with on each site a localized spin-^ 
state. A natural basis for the bilayer spin states consists of 
singlet and triplet states on each interlayer rung [12)[13]. 
The singlet state |0 0) equals ^ (| ti|2) - | Iit2)) and 
similar expressions hold for the triplet states |1 m) with 
m = -1,0, +1. 

In the Mott insulator the magnetism is fully determined 
by the interplay between the interlayer Heisenberg en- 
ergy J± and the intralayer energy J. For low values of 
a = J±/zJ (= J±/4:J on a square lattice) the ground 
state is an antiferromagnet, whilst in the limit of infinite 
J± the electrons on each interlayer rung tend to form sin- 
glets. Consequently, the system undergoes a 0(3) univer- 
sality class quantum phase transition at some finite J± 
from an antiferromagnetic phase to a 'singlet' phase. The 
zero-temperature phase diagram of the bilayer Heisenberg 
model 12 14 is shown in Figure [l]3. 

Following the premise that the Coulomb interaction is 
large compared to the kinetic energy, we can write down 
a low-energy effective exciton model where the doubly oc- 
cupied and vacant sites are all bound into an interlayer 
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J = 0.125 eV 
t = 0.1 eV 
a = 0.04 



QP dispersion 

bandwidth: t^/J ~ 0.1 eV 




Fig. 3: Expected exciton spectral function for the c-axis charge- 
transfer exciton in YBCO bilayers. We used model parameters 
J = 0.125 eV, t = 0.1 eV and a = 0.04. The exciton quasi- 
particle peak has a dispersion with bandwidth t^/J, and the 
quasiparticle peak is the most pronounced at the line between 
(7r,0) and (0,7r). Following at a distance of 2;t(J/t)^/^, a sec- 
ondary peak develops as a sign of Ising confinement. 



exciton state \E) consisting of a doublon-holon pair on a 
rung. This exciton can hop around, a process in which 
the exciton swaps places with the neighboring spin state 
described by the Hamiltonian 



<ij> \ 



0)i{0 0|,- + ^|lm)i{lm|,- ) {Ei\ 



(2) 

The sum runs over all nearest neighbor sites ij and the en- 
ergy t is the exciton hopping energy. Following the same 
strong coupling perturbation theory that led to the intro- 
duction of the t — J model based on the Hubbard model, 
we can express the exciton hopping parameter t = tl/V 
in terms of the bare electron hopping energy and the 
exciton binding energy V . The total Hamiltonian Ht-j of 
this new 'exciton t — J model' consists of both the exciton 
hopping terms ^ and the Heisenberg terms ([T]). 

Computations. — Using the exciton t — J model we 
can compute the exciton propagator 



G{k,uj) = {MEk- 



4lV'o) 



(3) 



where |?/^o) is the ground state of the bilayer Heisenberg 
model, and e\, creates an exciton state with momen- 
tum k. Next we compute the exciton spectral function 
A{k^uj) = — [G(/c,co')], which is directly related to the 
c-axis optical absorption HIs] 



Im [e-^{q\uj)\ - (t^fA{k = 0,cj) 



(4) 



where t±_ is the interlayer electron hopping. We expect 
that the exciton spectrum should be detectable at an en- 
ergy V below the electron-hole continuum, where V is the 
binding energy of the exciton. Note that the exciton spec- 
tral function can also be measured using Electron Energy 
Loss Spectroscopy (EELS) 16 or Resonant Inelastic X- 
ray Scattering (RIXS) 17 probing the c-axis. 



As can be seen in Figure [TJd, a moving exciton causes 
frustration in the previously perfect Neel state. In or- 
der to deal with this frustration in a manner similar to 
the methods developed for the single layer t — J model 
[5]|6], we first need to construct the correct linear spin 
wave theory (LSW) for the bilayer Heisenberg model 18 . 
Following [12^ we define the sum and difference of the spin 
operators as s = Si + S2 and s = Si — S2 respectively. 
The Heisenberg Hamiltonian ([T]) now reads 



(5) 



Next we introduce a mean-field approximation to the 
ground state on each rung by 



=7/,cosx|0 0), -sinx|10)i. 



(6) 



that interpolates between the Neel state (x = 7r/4) and 
the singlet state (x = 0), and rji = (—1)* alternates on 
the two sublattices. The angle x needs to be fixed by 
minimizing the ground state energy for a given a. On 
each rung there are three excited spin states orthogonal 
to the ground state 



6j 



(r?isinx|0 0)i+cosx|10)i)(G|i, (7) 
|1 l>i(G|i, (8) 
\l-lh{G\i. (9) 



and consequently we can express the spin sum and dif- 
ference operators exactly in terms of the spin excitation 
operators and 6^, 

5+ = V2r]i[cosx{bl^ - bi-) + sinx(^l+ei - 

si = -V^[sinx(^l+ + bi-) - cosx(^l+ei + ^i-ej)], 

sin2x(l-2eje,- ^6^6,,) 



-cos2x(e] + Ci) 
blb,^-blb,-. 



(10) 



We rewrite ([5| in terms of the spin excitation operators ^ 
[9]). The linear spin wave approximation now amounts to 
normal ordering the Hamiltonian and neglecting all inter- 



action terms between the spin excitation operators 19 



After performing a Fourier and Bogolyubov transforma- 
tion, we find the spin wave dispersions 



Jz 



(sin^ 2x + 2a cos^ X T 7k cos 2x)^ — 7^ (H) 



JzVC(CT7kCos2 2x) 



(12) 



with 7k = ^{coskx + cos ky) and C = sin^ 2x + acos2x. 
These dispersions correspond to the known spin wave spec- 
trum p!2|jT4] consisting of two single layer transversal 
(acoustic) spin waves (11) derived from the 6^ operators 
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Fig. 4: Feynman diagram representation of the Self-Consistent 
Born Approximation (SCBA). The self-energy of the exciton 
depends self-consistently on 'rainbow' diagrams where it emits 
and absorbs either one or two spin waves. The left two dia- 
grams contain interaction with the longitudinal spin wave, the 
right diagram contains the interaction with the transversal spin 
waves. Vertex corrections are neglected in the SCBA. 



and a gapped longitudinal (optical) spin wave (12) from 
e^. In the Neel phase, the transversal modes are gapless 
Goldstone modes while in the singlet phase all spin waves 
are degenerate gapped triplet excitations. 

Using the expressions for spin excitations ([7|-([9| we find 
that the exciton hopping term Q takes the form 

(ij) 

+ sin2x(el + e,) - ^6^6,,] + h.c. (13) 



Using the explicit transformations obtained by the Bo- 
golyubov transformation between local spin operators 
([9| and the spin waves we can rewrite the exciton hopping 
Hamiltonian into an exciton-spin wave interaction term. 
Subsequently we used the SCBA, which was successfully 
applied to the single hole problem [5j[6], to compute the 
exciton self-energy. In the SCBA one neglects the ver- 
tex corrections of the self-energy, which turned out to be 
a good approximation in comparison with exact diago- 
nalization studies [?]. Next to the fact that we need to 
include both types of spin waves, an extra complication 
arises since the self-energy also contains diagrams with 
two spin waves. The resulting diagrammatic expression 
for the exciton self-energy T^^k^uj) in the SCBA is shown 
in Figure |4] As for the single layer t — J model, the SCBA 
can only be solved numerically, and we have obtained the 
exciton spectral function A{k^uj) using an iterative pro- 
cedure with Monte Carlo integration over the spin wave 
momenta discretized on a 32 x 32 momentum grid. We 
start with self-energy S = and after approximately 20 
iterations the spectral function converged. 

Results. — The extensive details of the computations 
will be published elsewhere 20 and we will illustrate here 
the resulting exciton spectral function in representative 
limits. The most important control parameter is a since 
it is tuning the nature of the spin system, while t/J is just 
governing the degree of adiabaticity of the exciton motion. 
When a ^ 1, the magnetic ground state is just a stack of 
rung singlets which is effectively transparent for the exci- 
ton that propagates freely in this regime with a bandwidth 
of order of 8t. Given the expectation raised by the car- 
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Fig. 5: Exciton spectral functions when the interlayer magnetic 
coupling a = 0. In the top figure (t <^ J) the adiabatic motion 
of the exciton allows for the formation of a clear quasiparticle 
peak. In the anti-adiabatic limit (t ^ J, bottom figure) the 
spectral weight is shifted into a many-spin wave incoherent 
bump. Note that both spectra are almost non-dispersive. 



toon in Fig. [T]^, the most interesting limit is a ^ 0: here 
the antiferromagnetic order is strongest, while the inter- 
layer exchange coupling becomes inoperative with regard 
to flipping back [sj^ the two overturned spins caused by 
the hop of the exciton. In the strongly adiabatic regime 
where t <C J the exciton just localizes since it has to climb 
a very steep 'exchange hill' in order to move at all. Al- 
though its spectrum does contain all the 'ladder states' 
associated with its 'magnetic string' confinement these ac- 
quire a very small spectral weight because of the adiabatic 
limit: the first excited ladder state at energy Jz is just vis- 
ible with a spectral weight 0((fi / J'^). 

The confinement effect in this a ^ regime becomes 
fully visible in the anti-adiabatic t ^ J regime. Now the 
hopping exciton can explore the confining potential over 
large distances and the effect is that the spectral weight 
shifts to the highly excited ladder states leaving behind 
a strongly suppressed quasiparticle peak at the bottom 
of the spectrum. Both the adiabatic and anti-adiabatic 
results at exactly a = are shown in Figure [5| 

Whereas the ladder spectrum at a = seems to be ob- 
scured by a broad incoherent bump, inclusion of a finite 
a seems to intensify the ladder spectrum as is shown in 
Figure [2] To first order in a, the interaction of excitons 
with transversal spin waves is reduced. This favors the 
formation of a ladder spectrum, which can be seen if one 
neglects the transversal spin waves in the SCBA. The re- 
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maining interaction with the longitudinal spin waves leads 
up to first order in a to the following self-consistent equa- 
tion for the exciton self-energy, 



Jz — S(co' — Jz) 



(14) 



This is exactly the self-energy associated with a single hole 
in an Ising magnet [g] , where the naive picture of the spin 
frustration as drawn in Figure [T]^ becomes exact. Any 
motion of the dopant (electron, hole or exciton) away from 
its reference position increases the energy required for the 
spin mismatch. Oscillations around the reference position 
will generate a ladder spectrum, with the energy distance 
between the two lowest states scaling as t{J/t)'^^^. 

As compared to the Ising spin system the sharp ladder 
spectrum gets blurred due to the presence of transversal 
spin waves. The degree to which such blurring occurs de- 
pends on the ratio t/J. Consequently the ladder spectrum 
is most visible when t <^ J and at small but finite a, as 
can be seen in Figure [2] 

The ladder spectrum is, for example, expected for the 
half-filled cuprate bilayers in YBCO or Bi2Sr2CaCu208. 
Following earlier neutron scattering experiments on 
YBCO [l]|2T] one can deduce that the effective J = 125 ±5 
meV and J_\_ = 11 ±2 meV, which corresponds to an ef- 
fective value of a = 0.04ac where ac is the critical value 
of a 



12 



The question remains what a realistic estimate 
of the exciton binding energy is. The planar excitons are 
known to be strongly bound 10 with binding energy of 



the order of 1-2 eV. Since the Coulomb repulsion scales as 
V ~ (er)~^, we can relate the binding energy of the inter- 
layer excitons to that of the planar excitons. The distance 
between the layers is about twice the in-plane distance 
between nearest neighbor copper and oxygen atoms, but 
simultaneously we expect the dielectric constant Cc along 
the c-axis to be smaller than Cab due to the anisotropy 
in the screening. Combining these two effects, we con- 
sider it a reasonable assumption that the interlayer exci- 
ton binding energy is comparable to the in-plane binding 
energy. The hopping energy for electrons is approximately 
= 0.4 eV which yields, together with a Coulomb repul- 
sion estimate of F ~ 1.5 eV, an effective exciton hopping 
energy of t ~ 0.1 eV. The spectral function corresponding 
to this parameters is shown in Figure [3] Since t ^ J the 
ladder spectrum is strongly suppressed compared to the 
aforementioned anti-adiabatic limit. However, the Ising 
confinement still shows its signature in a small 'second 
ladder peak' at 0.4 eV energy above the exciton quasi- 
particle peak. To the best of our knowledge and to our 
surprise, the c-axis optical conductivity of YBCO has not 
been measured before. Detection of this second ladder 
peak in future experiments would suggest that indeed the 
interlayer excitons in cuprates are frustrated by the spin 
texture. 

Conclusion. — In conclusion, we have shown that 
the dynamical spin-hole frustration effects that are well 



known in the context of doped Mott insulators occur in a 
strongly amplified form dealing with interlayer excitons in 
Mott-insulating bilayer systems. We anticipate that the 
real significance of this discovery lies in the potential of 
creating exciton systems of this kind at a finite density. 
As compared to the usual doped Mott-insulators, excitons 
simplify matters because they are bosons rendering the 
situation much more tractable. At the same time, the 
tendency to exhibit complex ordering phenomena like the 
stripes and so forth are in first instance caused by dynam- 
ical frustration effects of the kind discussed in this paper. 
At least when the magnetic order is maintained, it ap- 
pears that the problem of excitons at finite density is sign 
free and can therefore be addressed in terms of bosonic 
field theory. We expect a coexistence of magnetic order 
and the exciton Bose condensate, but in the light of the 
present single exciton results it is clear that these will show 
unusual strong coupling features. 

This gives further impetus to the pursuit to create such 
finite density correlated exciton systems in the labora- 
tory. One can wonder whether such physics is already 
at work in the four-layer material Ba2Ca3Cu408F2 where 
self-doping effects occur creating simultaneously p and n- 
doped layers Much effort has been devoted to cre- 

ate equilibrium finite exciton densities using conventional 
semiconductors 23 , while exciton condensation has been 



demonstrated in coupled semiconductor 2DEGs 24 , 25 



In strongly correlated heterostructures, however, forma- 
tion of finite exciton densities is still far from achieved, 
although recent developments on oxide interfaces indi- 
cate exciting potential (see for example Besides 
the closely coupled p- and n-doped conducting interface- 
layers in these SrTiOa-LaAlOa-SrTiOa heterostructures, 
further candidates would be closely coupled p- and n- 
doped cuprates, such as YBa2Cu307_a; or La2-a^Sra;Cu04 
with Nd2-a:Cea^Cu04. The feasibility of this has already 
been experimentally demonstrated, e.g. in 27 , but the 
exact interface effects need to be investigated in more de- 
tail, both experimentally as well as theoretically ^28j|29j. 
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